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1 fH 1 ! Abstract 

In an attempt to solve Maxwell's first order system of equations, starting from a given 

c/j ■ initial state, it is found that a consistent solution depending on the temporal evolution 

of the sources cannot be calculated. The well known retarded solutions of the second 



O 



order equations, which are based on the introduction of potentials, turn out to be in 

^ \ disagreement with a direct solution of the first order system. 

O ' 
•g ■ PACS number: 03.50.De 

o • 1. Introduction 

In recent papers [1, 2] it was shown that Maxwell's equations have different formal 
solutions depending on the chosen gauge. In [2] it was argued that the formalism of 
gauge invariance is based on the tacit assumption of Maxwell's equations having unique 
£D ' solutions which appeared, however, not to be guaranteed a priori. In response to the 

publication of [2] it was pointed out in private communications [3] that uniqueness is 
a necessary consequence of the linear structure of the equations. These arguments are 
valid. If one finds, nevertheless, different solutions in Lorenz and in Coulomb gauge, it 
seems to indicate that a solution does not exist at all. Indeed, it was shown in [2] that 
the Lienard-Wiechert fields based on the Lorenz gauge do not satisfy the equations in 
the source region, unless one postulates a velocity dependent "deformation" of point 

J>V charges as in [1]. Furthermore, the formal solution for the vector potential in Coulomb 

gauge led to an undefined conditionally convergent integral which would even diverge 
upon differentiation. 

The reason for the difficulties encountered could have to do with the assumption of 
point sources which were exclusively considered in [2]. Therefore, it appears worthwhile 
to investigate the problem further, assuming smooth charge and current distributions 
as originally considered by Maxwell. In order to avoid any ambiguities arising from 
the introduction of potentials, it seems advisable to analyse directly the solvability 
of the first order system of Maxwell's equations (Sect. 2.). It turns out that the 
coupled first order system contains certain inconsistencies which prevent its solution 
when calculated by a numerical forward method proceeding in time. 

The usual method of solution derives inhomogeneous wave equations from the 
first order system, and expresses the solutions as retarded integrals by application of 
Duhamel's principle. In [2] it was argued that this method is not plausible, since the 
wave equations obtained by differentiating the first order system connect the travelling 
fields with the stationary sources at the same time, while in the retarded solutions the 
differentiation of the sources is inconsistently dated back to an earlier time. In Sect. 3. 
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we analyze the retarded solutions for smooth source distributions and find that these 
solutions do not satisfy the first order system. This is demonstrated in Sect. 4. by 
considering a specific example. 

2. The first order equations 

In vacuo the first order system as devised by Hertz on the basis of Maxwell's equations 
is supposed to describe the electromagnetic field: 

div Eg = 4n p (1) 

^ IdB 

rotE r = —— 2 

c at 

divB = (3) 

_ 4tt - ld(E g + E r ) 

rotB=—j + ^— >- 4 

c c at 

Here we have indicated that the electric field has two contributions of different struc- 
ture. In (^) only the irrotational part enters, whereas (J2J) contains exclusively the 
rotational part of the field. Both parts enter equation (0J). One may separate out the 
instantaneous contribution of the magnetic field and write (J3J) as two equations: 

4tt - 1 dE„ 
rotB = —j + -—f- 5 

c c at 

rotB 1 = — — (6) 

c at 

The quasi-static solutions of (JTJ) and (J3J) - subject to the boundary condition that the 
fields vanish at infinity - are represented by integrals over all space: 

E g (x, t) = ((J p (f , t) (x - x') -J™ (7) 

J J J \x — X \ 

x ,-, •> 1 [[f(-t,^i s 1 dE (x',t)\ ,_ _. d 3 x' , , 

B„ (f , t) = - /// (j (.?, t) + Tn ^j^j x (f - g) w - ¥f (8) 

It remains then to determine the rotational part of the electric field and the contribu- 
tion B\. 

Applying a numerical forward method one obtains from (jUJ) the difference equation: 

E r (At) = E r (0) + At c rot B x (0) (9) 

and from (J2J): 

B(At) = B(0)-AtcrotE r (0) (10) 

Assuming that the sources were constant for t < one has the initial conditions: 

E r (0) = By (0) = (11) 

Substituting this into (0) and (fTUj) one finds the curious result that neither E r nor the 
total magnetic field B proceed after the first time step, and this will remain so forever, 



at least in the vacuum region outside the sources. If the current would linearly rise to 
a new stationary level, e.g., equation (fTUj) would predict that B stays constant at its 
initial value, in contrast to (jHJ) which predicts that Bq rises simultaneously with the 
current and reaches a new stationary value as well. 
One may also split (J2j) into two equations: 

rotE rQ = - l - d A (12) 

c ot 

rotE rl = --^ (13) 

c ot 

The quasi-static solution of (fH?|) is: 

- 1 fff 8B x-x' 3/ , 1A . 

and from ()13|) follows: 

B x (At) = B x (0) - At c rot E rl (0) =0 (15) 

If B\ vanishes after the first time step as follows from (|15j) . and B stays also constant 
according to (fTUj) . a clear contradiction with (jHJ) arises. Furthermore, equation (0) 
predicts that the total rotational electric field stays constant, whereas the quasi-static 
part (fTlj) follows instantaneously all changes of B (t) according to (fPfj) . 

We note that the quasi-static expressions (J7|), (JBJ), (fTlj) can be seen as solutions 
of elliptic equations. On the other hand, one obtains from (J2J) and (@J) by mutual 
elimination of the fields the inhomogeneous hyperbolic equations: 



M-4^ = "-^+3^? (17) 



1 d 2 E r 4tt dj 1 d 2 E g 

c 2 dt 2 c dt c 2 dt 2 

As indicated in [2], the mixture of elliptic and hyperbolic equations inherent to Maxwell's 
system leads apparently to the inconsistencies which manifest themselves in the in- 
congruities implied in (JTHJ) as compared to (jBJ), and in ([1^1 as compared to ©• The 
system (1-4) does not permit a continuous temporal evolution from a given realistic 
initial state. In a region where the sources in (J16)) and (J 17)) vanish the homogeneous 
hyperbolic equations describe correctly propagating electromagnetic fields, but their 
production mechanism in connection to the sources remains obscure. 

Since in all textbooks it is claimed that Maxwell's equations do have solutions which 
are uniquely determined when the behaviour of the sources is given as a function of 
space and time, we must discuss the usual procedure to obtain these solutions which 
- according to our analysis - cannot satisfy the first order system. 

3. The retarded solutions 

The normal method of solution expresses the fields by potentials: 

^ 1 dA 

B = rotA, E = -V(f) — (18) 

Y c dt K J 



which leads to inhomogeneous wave equations in Lorenz gauge: 

1 d 2 6 

<»-*w = - 4 "' (19) 

A t 1 d 2 A 4tt -< 

**-*W = -T> (20) 

They are solved by application of Duhamel's principle to yield the retarded solutions, 
e.g.: 

A(x, t) = - [f[j&, t ~ & ~*\/c) rS^r (21) 

c J J J \x — x'\ 

Instead of introducing potentials one may solve the wave equations for the fields 
directly. The magnetic field, for example, can be expressed as the sum: 

B = B + B x (22) 

where B is the instantaneous part (JBJ), and B\ satisfies according to (fTSJl the equation: 

1 d 2 B x 1 d 2 B 

ABl " 7 2 ^v - 7 2 ~w (23) 

In analogy to (j2*Tj) this equation has the retarded solution: 

t> = t-\x -x*\/c (24) 



R, (t t) — fff 1 


(d 2 B (x',t>)\ 
y dt' 2 ) 


d 3 x' 


Bi{X ' t} 4k<?JJJ \ 


1 ~* ~*f \ 


Similarly, one may write: 







E = E + E 1 (25) 

and obtain from (|17|) a second order differential equation for E\\ 

. -. 1 9 2 ii 1 d 2 E , , 

where E is the instantaneous part of the electric field resulting from (jSJ) and ()12j): 

4w«-4///(^)w^r (27) 

The retarded solution of (|26j) is then: 

1 fff (d 2 E (£', t') \ dV 



It turns out that the fields as obtained from (|2*T|) and (|18|) are not the same fields 
as that calculated from ((22), ©, and (J2U), and from (J23J), (JUJ), and fl2EJ). This will 
be demonstrated in the next Section by choosing a specific example. Hence, we must 
conclude that the retarded solutions cannot be considered as true solutions of the first 
order equations. 

The reason for this failure must be sought in the inconsistency which lies in the 
fact that equations (|2L))1. ([23)1. (J2lj|) connect the sources j(x,t), B (x,t), E (x,t), 



respectively, with the travelling wave fields A (x, t), Bi (x, t), E\ (x, t) at the same time 
t, whereas in the retarded solutions (J2TJ), (J2U), (J25J) the differentiation of the source is 
dated back to the earlier time t' = t—\x — x'\/c. As pointed out in [2], the source may 
be very far away from the observation point, and may not even exist anymore when the 
fields B\ (x,t), Ei (x,t) are measured at time t. It makes little sense to differentiate 
non-existent instantaneous fields at time t, but this was necessary to derive equations 
dEj), (fTTjl from the system (1-4). Obviously, it constitutes a contradictio in adjecto 
connecting the travelling fields predicted by (|16|) and (|17|) with the stationary sources 
in the first order system at the same time. 

4. A specific example 

In order to facilitate the calculations we choose an example where we have div j = 0. 
In this case the scalar potential vanishes because of p = which makes Lorenz and 
Coulomb gauge identical: div A = 0. The chosen example is a hollow cylinder which 
carries a closed oscillating current driven by an rf-generator through a resistor R, as 
sketched in Fig. 1. It is assumed that the current was switched on at time t = — oo 
and oscillates with a sinusoidal time dependence: / exp (—iut). The current flows in 
a thin central filament, and returns symmetrically on the cylindrical surface. This can 
be achieved to an arbitrary degree of accuracy, if the inverse wave vector k = u/c is 
large compared to the dimensions of the device. 
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Figure 1: Oscillating current flowing in a closed circuit of cylindrical geometry 



The instantaneous magnetic field component (jHJ) produced in this configuration is: 
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— exp (— iut) , r < a , —b<z<b 
cr 



B 

B = 0, r > a , z < —b , z>6 
and the instantaneous electric field ()27|) becomes: 



(29) 



277 



cos (p' dr' 
R 



z'=+b 



z'=-b 



<,r dz^' = ° 



R 



ik I 
E = - — exp(-iut) I < 

2"7TC 



R = \jr 2 + r' 2 — 2 r r' cos tp' + (z — z') 
The retarded solution of the vector potential as obtained from (|2~Tj) is: 



> dp' (30) 



A 



Ie 



-iuit 



2-k I a 



2-kc 



I 



exp (i k R) 



R 



z'=+b 



4-6 



cos <£>' dr' e r + 



z'=-b 



exp (i k R) 
R 



r'=0 



cfe' e z > d<^' 



(31) 

It may be substituted into (J18JI to yield the fields as given by Jackson for a localized 
oscillating source [4]: 

27T ( a r „ • i, D ,„ . , „s , , ,n s=2+6 



5 



1 — ik R) s cos </?'dr' 




V (32) 



2-6 



r'=0 



£ 



ik I e 



-iujt 



2ir I a 



2ttc 



I 



e cos (p dr' 
R 



z'=+b 



e r + 



6 r „• 



z'=-b 



ikR dz l 



-b L 



/? 



r'=0 



> <V (33) 



where s = z — z'. It is doubtful whether these solutions satisfy also the differential 
equations PH|) and ()2b|). In order to check on this we consider, e.g., equation (J23j) 
adapted to our case: 



^.^-B^-AV,.^- *t>Jr ..,.., 



dr 2 dr x ^ ' dz 2 c 

where the right-hand-side must be set to zero outside the cylinder of Fig. 1. We 
integrate this equation with respect to r and obtain: 



dB x B x 1_ 



'd 2 B 1 

dz 2 



k 2 B x dr 



k 2 I 



-ILUt 



(35) 



The contribution B\ may be calculated from (32) by expansion of the exponential 
function for kR < 1. In zero order one obtains the instantaneous field (|29|h and in 
second order one has: 



Bi 



Ik 2 e 



2„—iut 



L'TT 



An c 



s cos ip'dr' 
R 



= 2 + 6 z + 6 



= 2-6 



(r — r' cos <p') dz' 
R 



r'=0 



<V+0 (fc n>2 ) 
(36) 



The integration over r' and z' may be carried out analytically to yield: 



Ik 2 e~ iuJt 2t 



1 fee r i— 

Bi = — / [s cos ip' In (r' - r cos </?' + R) + (r — r' cos <p') In (s + R)] s s Z z z _ b yZl dip' 

o 

(37) 

Expanding this expression in a power series of r, and inserting it into the left-hand-side 
of flMJ) we find for z = 0: 



JPe 



-iidt 



2c 



(a 2 + 6 2 ) 




+ O (r n>4 ) (38) 



which is obviously at variance with the right-hand-side of ()35j) . A similar conclusion 
is reached, if (155)1 is substituted into (|26|) . This can only be checked numerically, since 
the instantaneous field E does not vanish outside the cylinder, in contrast to B . 

Result (|3*8*j) proves that the standard solutions (}3*2"j) and (|3*3*j) do not satisfy the first 
order system from which equations (fTBj) and (fT7|) were derived. Hence, our conclusion 
in Sect. 2., namely that the first order system does not permit a solution, cannot be 
refuted by referring to the retarded solutions as taught in the textbooks such as [4]. 

There is also a physical reason to reject Jackson's solution (|31|) for the considered 
case. If one calculates the fields with ()18j) from (j31|) and evaluates the Poynting vector 
E x B at large distance, one can integrate the total radiation power emitted by the 
closed circuit of Fig. 1: 

«.-U±.W.*-M* ,30, 

This result is obviously not physical. The device in question may be seen as a short- 
circuited cable which should not continuously loose energy to the outside world; in 
particular not when the enclosing shell would be made out of superconducting material. 
The predicted power loss (J3l?j) could certainly not be confirmed experimentally. 

5. Conclusions 

It has been shown that an attempt to calculate numerically the temporal evolution 
of the electromagnetic field from the full set of Maxwell's first order equations will 
fail due to the internal inconsistencies built into the coupled system of equations. As 
noted earlier [2] , the reason lies in the fact that the travelling wave fields are connected 
with the stationary sources at the same time. 

Maxwell's equations describe correctly the production of the instantaneous electro- 
magnetic field, and also the propagation of wave fields in empty space. The production 
mechanism of electromagnetic waves by time varying sources, however, does not find 
an explanation in the framework of Maxwell's theory. Contrary to what is commonly 
believed, the retarded solutions for the electromagnetic potentials do not lead to fields 
which are in agreement with a direct solution of the second order differential equations 
for the fields. 
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